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Abstract. We obtain an asymptotic solution for £ of the Cauchy prob- 
lem for linear first-order symmetric hyperbolic systems with oscillatory initial 
values written in the eikonal form of geometric optics with frequency 1 /e , but 
with complex phases. For the most common linear wave propagation models, 
this kind on Cauchy problems are well-known in the applied literature and 
their asymptotic theory, referred to as complex geometric optics, is attracting 
interest for applications. In this work, which is the first of a series of papers 
dedicated to complex geometric optics for nonlinear symmetric hyperbolic sys- 
tems, we develop a rigorous linear theory and set the basis for the subsequent 
nonlinear analysis. 



1. Introduction 

We shall obtain an asymptotic solution for e — > of the Cauchy problem for 
linear first-order symmetric hyperbolic systems in several spatial dimensions with 
initial values in the form, 

m 

h^{x) — ^ hf^{x)e^^''^^^^^ , (V'm(^) complex-valued phases), 

that is, a superposition of waves of frequency 1/e, amplitude and complex phase 
ijj^ for which we always assume the condition Im^ijj^) > 0. 

More specifically, we construct a family of waves such that — t^j^j^gl ~ O(e^) 

and \L{t,x,d)v^\ = 0{e^) where L{t,x,d) is the considered symmetric hyperbolic 
system; then, we prove the convergence of the asymptotic solution to the exact 
one in a suitable topology. 

The strategy will be introducing complex-valued phase functions in the frame- 
work of geometric optics for which we refer to the lectures by Ranch ^ where 
rigorous results are given together with a comprehensive list of references. In the 
applied literature [2-6] the theory of oscillatory waves with a complex-valued phase 
is widely developed in several different variants, cf., the recent book by Kravtsov 
for a tentative classification; particularly, the beam tracing method of Pereverzev 
|3[ ^ is currently used in fusion experiments and complex geometric optics finds 
several applications in electrodynamics and geophysics |^ ^. On the other hand, 
to our knowledge, a rigorous justification of such approaches is still lacking even in 
the linear case, especially for systems, and no attempts have been made to address 
the nonlinear case. Here, we put the complex geometric optics theory on a rigorous 
basis by means of novel tools that clarify the analysis of the polarization in the 
presence of complex phases as well as by a precise argument based on a partition 
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of unity. Although it may be of independent interest, this result is meant to be 
the first part of a work dedicated to the study of such initial value problems for 
nonlinear symmetric hyperbolic systems. 

The main point is that the non-negative imaginary part of the phase yields a 
strong localization of h'^ around the set R° = [J^{x; ImV';^(a;) = 0}, as e ^ 0, thus, 



h'^ exhibits features common to both wave trains (oscillatory character) and wave 
pulses (localization): the more the wave field is localized, the more it oscillates in 
such a way that, qualitatively speaking, the number of oscillations under the wave 
field envelope can be thought of to be independent on e. For such localized waves, 
diffraction effects should be taken into account in the leading order asymptotics, 
even for short-time propagation. Our analysis parallels the results on short pulses 
(?[ ^, albeit diffraction is accounted for in a different way, i.e., through the coupling 
of the real and imaginary parts of the phases. Such a description of diffraction 
should be compared to a deformation of the wave field envelope described by a 
Schrodinger-type equation relevant to the paraxial approximation (also known as 
"parabolic wave equation" |^). In the forthcoming papers we shall see that this 
duality between complex phases and "Schrodinger-type" evolution of the wave field 
envelope is still valid in the nonlinear context. 

2. General Assumptions and Main Results 
We shall consider systems of first-order partial differential equations of the form 

d 

L(t,x,d)u{t,x) = dtu{t,x) + ^A^(t,x)dx^u(t,x) + B{t,x)u{t,x) =0, (2.1) 

for the wave field u € C°^{il;<C^) on the closure fl of an open, connected and 
bounded set n C IRi+'*; these are the systems considered by Lax in his seminal 
paper [|o). Here, Aj, B are functions valued in the space End(C^) oi N x N 
complex matrices depending smoothly on {t, x). 

Assumption 1. The matrices Aj{t, x) are Hermitian in a neighbourhood ofil. For 
the particular case of a single equation (N — \), this reduces to Aj{t, x) G R. 

Assumption |l| is equivalent to requiring that the matrix 

d 

A(t,x,C) =^A,(t,a;)e,, 
i=i 

is Hermitian for {t, x) in a neig hbourhood of Q. and fe R'', with R'' = R'' \ {0}. 



This implies that the principal part of the operator ( 2T ) , namely, 

d 

L^{t,x,d)=dt+Y,Aj{t,x)d,^, (2.2) 

is strictly hyperbolic in time, that is, the characteristic equation detuLg ~ 0, with 

aLo {t, x,T,^) = i{Tl + A{t, x,£,)), 
being the principal symbol of Lq, has real- valued roots given by the solution of 

fl{t,X,T,^)=T + Xl{t,X,O=0, 
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with A; the l-ih eigenvalue of A. This defines a conic variety Char(L) in T*R^+'^ \ 
0, the cotangent bundle without the zero section, which is called characteristic 
variety or local dispersion surface in the physics literature; we will assume that 
different branches of the characteristic variety corresponding to different eigenvalues 
A; cannot merge; specifically, we have the following assumption. 

Assumption 2. The eigenvalues Xi{t,x,^) have constant multiplicity. 

We know that eigenvalues of constant multiplicity are everywhere distinct, thus, 
one can label them so that Ai < A2 < • • • < A/ < • • • , and, through the implicit 
function theorem, one can show that they depend smoothly on the entries of the 
matrix A [|ll]; in addition, they are homogeneous of degree one, as follows by 
inspection of the characteristic equation. Moreover, there is a constant r > such 
that \Xi{t,x,^) — Xi'{t,x,£^)\ > 2r for {t,x,£,) e 17 x S'^~^, hence, we can use the 
contour integral representation for the spectral projectors, [Q, |l2[ , 

Mt,x,0^:^ <f _ {zI~A{t,x,i)y^dz, (2.3) 



which shows that TTi{t, x, ^) = TTi{t, x, ^) are of class C°°(S1 x M,'^) and homogeneous 
of degree in the variable ^ = |^|^. Finally we give the definition of complex phase. 

Definition 2.1 (Complex phase). Form > integer, (j> € C°°{0; C™) in a bounded 
open set O C R" is a complex phase iff, for each component 4>^ = </5^ + X/x; one 
has dif^i =/= and > m C 

Now we choose the domain fl according to 

n = {{t,x) e R^+'^; <t <T,\x - x\ < p- ct}, 

for suitable constants x £ M!^, p > 0, T G (0, p/c), with c > being the propagation 
speed of the operator {2A); we recall that such an f2 is a domain of determinacy 
for the operator L and the propagation speed c is characterized by the condition, 
0, 

d 

j=l 

on the boundary of Ti; we denote X*' =Un {t = t'}, for t' £ [0,T], with X* the 
interior of X*. Then, we consider the Cauchy problem for equation ( ^.l[ ) with initial 
value 

ni 

4^o(a;) = h'ix) = h^i^V''-^"^^', ^ e X^, (2.4) 

where e G 1R+, the functions ip = (ipf^) G C°°{X°;<C"^) are complex phases and 
hf^ € C°°(X°;C^) are called amplitudes. We require that the initial complex 
phases fulfill the following geometric hypothesis. 

Assumption 3. The locus l'ni^f^(x) — amounts to a closed submanifold i?° C X° 
without boundary and i?° n R° ~ for p ^ v so that R° = IJ^ i?° is also a closed 
submanifold without boundary. 

As for the amplitudes /i^ G C°°(X°;C^), we assume the following condition 
which can always be satisfied after splitting each as appropriate. 
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Condition 1. For every ii there is I ~ l{fi) such that 

TTi{0, X, dReTl)f_,{x))h^{x) = h^{x), 

which means that each component of the initial datum is polarized in the eigenspace 
of a specific eigenvalue A/ . Such specially polarized waves are referred to as eigen- 
modes in the applied literature. 

We search for an asymptotic solution G C°°(ri; C^) in the form 

v'{t,x) = ^a^(^,x)e^^-(*^-)/^ al{t,x) = a|]^) (t, a;) + £a« (i, x), (2.5) 

with (j) — {4i^)fi being a multi- valued complex phase and a|f'', j = 0,1, are the 
amplitudes. Then, we readily get 

L{t,x,d)v^ = ^ |^e"VLo(t,a;, d(/)p)a^°^ + (jLo{i,x,d(j)^^)a''jl^ 

+ Loit, x, 9)a(f ) + Bit, x)a(f ) + eLoit, x, 9)aW + £B{t, x)a'j}A e'^-/^ (2.6) 

where the principal symbol is a polynomial in (r, ^) and it extends to an entire 
function in such variables, henc e, ULo{t,x, d4>^) makes sense for complex- valued. 

Usually, one imposes in (2_^) that the coefficients of different powers of e vanish 
identically. Here, instead, the presence of the exponential e~^i^/'^ , — Imt/)^, 
allows us to consider weaker conditions; precisely we have the following estimates 
used independently by Maslov and Pereverzev 

Lemma 2.2. Let k > be an integer, (1) = ip + ixe C°"(f7; C), / £ C°°(n) and let 
S d be any (non-empty) set such that S D {x{t, 2;) = 0} = 0. Then, 



£ 



-k 



^|/e**/^| < Ck, for every (t, x) e S and e eR+, 



where Ck = k^e ^ &n■p^^^^^)^-g\f{t,x)/x^{t,x)\. 

Proof. Since x > in S, the function e 1— > e^'^'e^-^^^ is bounded for every (t, x) ^ S 
with maximum value [k/ex)^ . D 

Hence, and Lv^ are localized around the set 

R=\JR^„ with i?^ = {(i,x) e fl;x^(t,a;) = Im0^(i,a;) = 0}; 

both R and i?^ are assumed to be submanifold of 17 characteristic for the operator 
( 2.1) and referred to as the reference manifold for the complex phases (j) and 0^, 
respectively. 

Assumption 4. In a conic neighbourhood of in T*]R^+'^ \ 0, there are isotropic 
suhmanifolds Ai, . . . , Am such that, for 11 £ {1, . . . , m}, 

(i) {0,x,T,£,) e A^ only if 

X e R°^, ^^dil;^{x), 

with i?° given in assumption^ and ip^ given in the initial value \2.^ ; when 
X S dtp^ix) is real valued; 
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(ii) with /(/i) given in condition^, and for every {t,x,T,£^) G A^, 

fi(M){t,x,T,£,) = 0, and Hf^^^^^ e T(^t.x.T.i)^f^, 

where Hf is the Hamiltonian field corresponding to the Hamiltonian func- 
tion f ; 

(iii) the canonical projection T*M.^~^'^ — > R^+'' restricts to diffeomorphisms be- 
tween n {(t, x, r, (t, x) G ri} and a closed submanifold i?^ C 51; 
dimi?^ — dimi?° + 1 and is transverse to ; 

(iv) the submanifolds i?i , . . . , Rm are disjoint. 

Remark 2.3. The construction of the isotropic submanifolds A^ is exactly the same 
as in the standard geometric optics method, the only difference being its dimension. 
With minor modification we can also consider the (degenerate) case in which i?° 
amounts to a single point and, correspondingly, i?^ is a single geometric optics ray; 
this is the case considered in most applications [3-6]. 

Assumption ^ addresses the two main geometric objects of the theory, that is, the 
isotropic manifold A^ and the corresponding projection i?^. In the theory of Maslov 
(D the focus is on the former, whereas in the paraxial approach by Pereverzev ^ 
one looks at the latter. In section |4[ we shall see that assumption ^ is, indeed, 

It allows us, in 



2.1 



an hypothesis on the characteristics flow for the operator 
particular, to construct a coordinate chart (O^, k^), where is a neighbourhood 
in O of a point (to,a;o) € i?^ and : 9 {t,x) ^ {t,r,s) £ [0, T] x x Os 
is a diffeomorphism, with Or C Os C 'Ef^ {di -\- d2 = d), so that n i?^ is 
mapped into {s — 0}; the latter is just the submanifold property and we see 
that the boundary 9-/?^, corresponding to t = and t = T, is plainly accounted for. 

In addition, in section |^, we shall see that XtJ.{t^^) ~ Ini(/)^(i,x) > c|s|'' for a 
constant c > and an even integer q > 0. 

Lemma 2.4. Let (j) and f be as in lemma and let R = {Imcj) = X 0} be 
a submanifold admitting coordinates {t, r, s) = K{t, x) on a neighbourhood O C f2 
as described above. We assume further that x(t,r,s) > c|s|' in [0,T] x Or Os 
where c > and q G ¥i is an even integer. Then, for every fc G M, t e [0, T] and 
A C Or X Os open and bounded with A n {s = 0} 0, there are constants Ca, 
a e INf''^ with length \a\ = k, such that 

\{f{t,r,s) - Ca(i,r),s")e^*(*''-'^)/^| <£^ ^ sn^Jd'^ f{t,r,s)\C^, 

\a\<k \a\=k {r,s)(^A 

uniformly for (r, s) G A, £ G R-j,; such an estimate can be made uniform in [0, T] x 
A. Here, flt,r,s) ^ f o K^^(t,r, s) and analogously for the other functions, whereas 
c^it,r)=dffit,r,0)/a\. 

Proof. First, let us fix a time t G [0, T]. By Taylor formula we have 
f{t,r,s)~ Ca{t,r)s°' ^ Y s"ca(t,r,s), 

|a|</i: |a|— /c 

where \ca(t,r,s)\ < sup^ r, s)| in A. Moreover, if |a| = fc. 
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with Cq being the maximum of |ti"|e^'^l"l'' for v e R''^. This yields the claimed 
estimate pointwise in [0, T] which, on the other hand, entails the miiform estimate 
in [0,T] X A e.s the right-hand side is continuous on [0,r]. □ 

In virtue of lemma ^.2| we see that it is enough to consider a neighbourhood of 



the reference manifold R and lemma 2.4, specialized for (j) = (p^^, tells us that it is 
enough to find a solution of 

ai„(t,a:,d0^)a(f) =0(|s|^), (2.7) 

aLo{t,x,dq^^)al'^ + Lo{t,x,d)al"^ + B{t,x)a^^^ ^ 0{\s\i), (2.8) 

in a neighbourhood of R^^ . The lower non-trivial value of the parameter q is clearly 
q = 2 and this is the case we are interested in. 

The first result of the paper is the existence of an equivalence class of solutions 



to equations (2/7) and ( |2.§| ). This can be done by means of standard tools: an 
extension argument put forward in section ^, the construction of an approximate 
solution for the so-called complex eikonal eq uati on for addressed in section ^ to- 



gether with the analysis of equation (2/7) and (2^) in section |5. The construction of 
the complex geometric optics solution is completed in section | where the following 
proposition and its corollary are proved. 

Proposition 2.5. Let assumptions [J-^ be satisfied together with condition [J and 
d1hmpfj_{x) > when x G Then, there exists an equivalence class of functions 
= J2 a^^e*"*"/^ e C°°(n; C^) such that, for e £ (0, Eq], < Eq < 1, 

a) j/i"^ — t^j^f^pl < CiE"^ , uniformly in X°; 

b) \L(t,x,d)v^\ < C2S^ , uniformly in Vt. 

For every fixed e G R+, one can apply the classical existence and uniqueness 
results based on the energy integral method ^ that gives as a byproduct the 
convergence of the complex geometric optics solution to the exact solution. 

Corollary 2.6. Let u^,v^ G C°°{fl) be the exact solution and any representative 



of the equivalence class in proposition (2.L), respectively, then 



sup \\u''[t) -v''{t)\\L2^x^) <Ce^ , /ore g(0,£o]- 

0<t<T 

In view of the specific construction of the approximate solution, which is based 
of the smooth reference manifold R, it is expected that the latter L^ estimate can 
be refined to an estimate by means of linear conormal estimates; for instance, 
the case of a single wave (m = 1) with a codimension one reference manifold is 
the analogous of the case considered by Alterman and Ranch in their study of 
nonlinear geometric optics for short pulses ||^. However, we delay the study of the 
appropriate conormal estimates to the work on generic nonlinear systems. 

3. Analysis of the Matrix aLo{t,x,d(j>) with a Complex Phase 

As a preliminary analysis we study the kernel and the range of the matrix 
aLo{t,x,d(j)) for (/> G C°^{il;C) being a single complex phase. 

Although we know everything about (TLoit,x,dip) for real-valued cp, the cor- 
responding results for a complex-valued phase 4> do not follow directly from the 
assumptions. On the other hand, we note that \dx\, X = Imc^, is small near the 
zero level set R — {{t,x) G r2;x(^,a;) — 0} since x restricted to any curve 7 C 
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has a minimum in 7 n i? when 7 is transversal to R and it is constant if 7 C i?; 
hence it is natural to regard idx as a perturbation in the expression dcf) — dip + idx- 
For any function / G C°°(R'*;E) taking values in a generic finite-dimensional 
vector space E and for any integer n G IN, let us set 

/'"^(C) = E '-^d^fiOv", C = e + e e ^ 0. (3.1) 

^ — ^ a! ^ 

I Q I < n 

We see that is of class C°° in the complex domain R'^ + iR'^ C C''. 

Proposition 3.1. Let e C°°(R'';E,), i = 1,2,3, and let us assume that there 
exists a composition law : Ei x E2 9 (/i,/2) ^ fs = /1/2 S E3 such that the 
Leibniz's rule d^-{fif2) = (9j,/i)/2 + fi{d^i.f2) holds; then, for every rt G IN, 

with p^'^\i,ri) = J2 ^^9|'/i(09f/2(0?7"+'' = 0(|77|"+i), the sum being over a, (3 
with n < |a + /3| < 2n. 

Proof. By applying the Leibniz's formula we have 

|7|<n |7|<''^ a+/3— 7 

whereas, 

/(")(c)i^")(c)^ E E '^d-ih{Odlf2{or^' 

|a|<n|^l<" 
, „-|a+/3| 

\a+l3\<n 

□ 

Remark 3.2. Clearly the definition of /'"^ and the corresponding property for 
composition laws still hold if f depends also on {t,x) G fi. 

We now apply this simple result to the matrices A{t,x,£^) and 7r;(t,x,^) defined 
in section |[ the composition law being the matrix multiplication. 

Proposition 3.3. For every n G IN, 

a) the exact identity / = X); ^["^ holds; 

b) the following identities hold modulo 0(|?7|"+"'^), 

c) there is a unitary matrix U(i,x,^) such that, modulo 0(|77|"^""'^), 

/ = u(")v(") = v(")u("\ u(»)I(»)v(") = diag(Aj"\ . . . , A^'), 

where V = U* and diag denotes the diagonal matrix, the eigenvalues a|"'* 
being counted with their multiplicity. 
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Proof. We know that / = tt/, tt/ttc = Swtti, A — A/tt; and Att/ — Xini, then 
we have also = SI'tt; for every multi-index a with \a\ > 1. The latter implies 



/ |Q|<n ■ I I 



which is the exact identity in a) . The relations in b) follow directly form proposition 
3.1, for instance the first one reads 

and analogously for the others. As for c), we know that A is Hermitian, hence, 
there exists a unitary matrix U(t, x, ^) such that UAU* = diag(Ai, . . . , Xn) where 
the eigenvalues are counted with their multiplicity. Then, the claimed identities 



follows from proposition 3.1. □ 
Corollary 3.4. For every integer n > and complex phase (f>, 

-taL„{t, X, d<j>) = (dtcl) + (t, X, 40))^=}"^ {t, X, d^cb) + 0(|d,.xr+^), 
U(") (-iai JV(") = diag(9t0 + A^") {t, x, d^^), . . . , S^c/. + A^^' {t, x, d,<l>)) 

+ 0{\d^xr+'). 

Proof. Since d^A{t,x,^) = for |a| > 1 we have I(")(t,a;,C) = A{t,x,C) and, by 
definition d(l){t,x) e + iM.'^; thus, -iaLoit, x, d(j)) = dt(t>I + A'''^\t,x,d^(t)) 
for every n>l. The claim then follows from proposition 3.3. □ 



Corollary 3.4 gives enough informations on the kernel and the range of the matrix 
aLo{t,x,d(j)) modulo the 0{\dxX\^~^^) remainder. In particular, we see that the 
kernel is non-trivial if we pick (j) such that 

(i, X, d^) - dt^ + a|"' (t, X, d^cj)), 

is 0{\dxx\^'^^) for some I, that is, if we find an approximate solution of the complex 
eikonal equation 

D^,^\t,x,d(l){t,x)) = 0, 
within a fixed order of accuracy. 

4. Approximate Solution of the Complex Eikonal Equation 

We will now determine the approximate solution (f) G C°"{Vl\ C) of the Cauchy 
problem 

J9t0-|-A(")(t,a;,d:,0) =0, (i, G H 

with X{t,x,£,) e C°°{n X R"^), homogeneous of degree 1 in ^, and 2p e C°°(X°; C). 
The analogous of assumptions || and ^ with rn = 1 and f{t, x,t,^) = t + A(t, x, ^) 
are supposed to be true. 

According to assumption ^ there is an isotropic submanifold A which is deter- 
mined by the fiow-out along the Hamiltonian vector field Hf in T*Ii^^'^ of the 
isotropic submanifold given over the hyperplane {t = 0} by the graph of dRe^p(x) 
for X E R°; the projection of A onto amounts to the submanifold R which is 
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supposed to be transverse to the boundary of the domain il. If the map- 
ping F : [0,T] X R" ^ il is the projection on il of the Hamiltonian flow, then 
R ~ F{[0, T] X R°) and we have the following proposition. 

Proposition 4.1. If assumption^ is satisfied, 

a) for every t £ [0,T], i?* = X* Cl R is a smooth submanifold of X* and 
i?* ^ R°; 

b) R is fihered over R° , the fibers being the integral curves of the smooth vector 
field V{t^x) = dt + d^X{t,x) restricted to R. 

Proof. According to assumption U, the Hamiltonian flow amounts to a diffeomor- 
phism : [0, T] X i?° ^ A n {{t, x, r, 0; {t, x) G H} and A n {(t, x, r, 0; (i, x) e Tl} is 
diffeomorphic to i?, thus, F is an embedding of [0, T] x R° into and F^^ is well- 
defined on i?; specifically, x) = (t, with x — x{t, x°) being the projection 
of the Hamiltonian orbit issuing from {x° ,dip{x^)) (d'ip{x) is real when x 6 R°). 

a) The mapping F* = F{t, ■) : R" A* given by x° ^ (t, x{t, x°)) = F{t, x°) is 
an embedding of R° in fl for every t e [0, T], therefore i?* = F*(_R°) is a smooth 
manifold without boundary diffeomorphic to R". 

b) On fixing some local coordinates r in an open neighbourhood TZ° C R° we 
have natural local coordinates in the neighbourhood TZ = F{[0, T] x 7?,°) C R given 
by {t,x) {t,r), with r defined as the coordinates in 72.° of the point x° such 
that X — x{t,x°). Moreover, we have the projection H : R 3 {t,x) ^ G R°, 
where, again, x and x° are related by x — x{t,x°); in coordinates H amounts to 
the projection on the second factor {t,r) i-^ r. The fibers of H are Il~^{x°) — 
{{t,x{t,x")) G R;t G [0, r]}, that is the integral curve of the vector field V. □ 

We shall now extend the fibered coordinates (t, r) to coordinates (i, r, s) in a 
neighbourhood O C fl. With this aim, we make use of the Euclidean metric (•, •) 
in order to define the normal space of i?* in A* at the point {t,x) G i?*, namely, 

iV(,^,)i?* = {Sx G T(,,,)A*; {dr,dx) ^ 0,V<5r G T^t,,)R'}, 

and the disjoint union for {t,x) G i?* of A^j-j aj-ji?* defines the normal bundle NR*. 
Then, the disjoint union 

T= [j TXIj^,, AA= y NR\ 

0<t<T 0<t<T 

are smooth vector bundles over R locally isomorphic to i? x H'^ and R x IR,''^^ 
respectively. To see this, let e° G N^oR°, 5r° G T^oR" and Sr = DF*Sr° G 
Tf-t(2.o-)i?* with DF* be the Jacobian matrix of the diffeomorphism F* : i?° — > i?* C 
A* for i G [0,T]; then, 

0= (e°,(5r°) = {e°,{DFY^Sr) = (*(i?F*)-ie°, (5r), 

and defines a local bundle isomorphism : NR° NR* for all t G [0,r]. 

Therefore, one can extend any vector bundle chart on A^i?° to Ai?* smoothly for all 
t G [0, T], thus, obtaining a vector bundle chart for of the form [0, T] x 72° x R''^. 
as for T, we can construct a vector bundle chart on noting that jj-jA'^t = 
T{t,x)R^ ® N (j^ j.)R^ and on applying the same argument to the relation 5r = DF^5r° . 
Furthermore, A/" is a subbundle of T and let F be the projector : T — > A/"; when, 
evaluated in a point (t, x) G i?, F(i, x) amounts the projector : T(^t,x)X^ A^(t,a;)-R*- 
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An orthogonal frame for M can be conveniently obtained by solving the Cauchy 
problem 



d2 

[TdtT + T)dt{I - r)]e, = c»jre,> 

j=i 



with orthogonal initial data, namely, (e°,e°) = 5ij, and with Cy S C°°(i?;R), 
+ = 0. Since T(t.x)X^ = K-'^, N(t.x)R^ C E'*, we have T{t,x) e End(R'*) and 
ei G R''; indeed, that amounts to a system of linear ordinary differential equations 
in M,'^ for which the solution exists in [0, T] and it is unique. On the other hand, 
we see that, if {ei}i is a solution, also {Tei}i does, hence, it should be = Fe^ 
or equivalently, e C°°{R\N) as required. Finally, such a solution defines an 
orthogonal frame for TV. The orthogonality can be proved on considering the matrix 
E = {Eij) whose elements are just the scalar products Eij = {ei,ej). If e; is a 
solution, then dtCi = [2F(9tF) - c^tFje^ + J^j^^v^^^j = [dt^i^^i + Y^j^^ij^j and, 
on noting that the commutator is anti-symmetric, i.e., *[9tF,F] = — [9tF,F] that 
follows from *F = F, one finds 

dtE = CE + E *C, £^14=0 = I, 

where C = (cy) is anti-symmetric. By direct substitution we see that E{t^ r) = / is 
the unique solution proving the orthogonality. One should note that the construc- 
tion of the basis {ei{t,r)} does not depend on the coordinates r on R°, i.e., the 
projector F(t, r) and vectors etit, r) are invariant under a coordinate change in R°. 
We can now define coordinates [0, T] x x Cs 3 {t, r, s) ^ {t, x) E il hy 

X = x{t, r, s) = x{t, r) + ej{t, r)sj, (4.1) 

for suitably small neig hbourhoods Or C R''i and Os C R'*^ We see that the 
differential of the map (M.ll) evaluated at s = amounts to 



which has a non-zero determinant as the columns are linearly independent. There- 



fore, the map (4.1) is a local diffeomorphism in a neighbourhood O, with TZ C O C 



fl and with boundary dO ^ 0° U O"^ , where 0* = O n X* are open in X*. These 



are the coordinates used in the formulation of lemma 2.4 



Remark 4.2. For the degenerate case R° — {x°}, the coordinates s are simply 
given by s — x — x{t,x°). However, one can also make use of the Frenet frame 
associated to the curve R, j^] . This differs from the construction of local coor- 
dinates described above in the fact that our normal vectors Ci lie in t — constant 
hyperplanes. 

We will make use of such coordinates to find an asymptotic solution of the 
Cauchy problem for the complex eikonal equation. Specifically, let us set 



|a|<n 
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the sum being over the multi- index a € IN''^; we want to choose the coefficients 
^a{t,r) so that 

• the manifold B = F{[Q,T] x R°) is the reference manifold for (p, 

• the complex eikonal equation is satisfied near R modulo 0{s^~^^), i.e, we 
have 

aj'£)(")(i,r,s,#(t,r,s))|^ = 0, for |a| < n. (4.2) 

Here, the subscript \R means that the function should be evaluated for s = 0. 

We see that (j)o{t,r) = (j)\i{{t,r) whereas for |a| = 1, amounts to (j)j{t,r) = 
9sj<f>\R{t,r) and, since R should be the reference manifold of </>, we must have 
(j)o{t,r) = ipo{t,r) e R, and cf)j{t,r) = ds,(p\nit,r) £ R, for j G {1, . . . ,^2}. Then, 
the lowest non-trivial order, for which complex-valued phases are found, is n = 2, 
namely, 

(j){t,r,s) = (poit,r) +^(pt{t,r)si + - ^ (j)^j{t,r)s^Sj, (4.3a) 

i= 

and we shall make use of notations 

Re(p{t,r,s) = (p{t,r,s) = (po{t,r) + '^(pi{t,r)si + ^ X] V>ij{t,r)siSj, (4.3b) 

d2 

2 



d2 



i=l 



^ d2 

lm(l){t, r, s) = x{t, r,s) = - ^ Xij{t, r)siSj. 



(4.3c) 



The matrix {xij ) should be positive definite. As a direct consequence we have. 



drk 



dfo d(j) 



R drk ' dsj 



R — drkdrh 



R drkdrh ' 



drkdsj R drk ' dsidsj 

From now on we consider only the lowest non-trivial order and we drop the index 
n = 2. First, we note that, for every t £ [0, T], the function 

-ViVj 



~Xit, X, C) = A(, X, + . E -If'^'' ^' ^) ^ 



dm^ 



is invariantly defined if {x, ^, rj) G T*X* ©T*X*, where © denotes the Whitney sum 
of vector bundles. Thus, we can readily write X{t,x,dx(l>) in terms of coordinates 
{t,r,s), namely. 



X{t,x,d(l){t,x)) = A{t,r,s,dr(p,dsip) +« V^t;— t; — (t,r, s,dr(p,ds(p) 

^ drk dpk 



dx dA 
dx dx d^A 



E 



dsi da 
dx dx d^A 



2 ^ drk drh dphdpk 



(t, r, s, drip, dgip) 



■{t,r,s,dr(p,dsip) - ^E 



dx dx, d^A 



drk dsj dpkdcTn 2 ^ dsi dsj daida-j 

kj J > J J J 

with (p, a) are the coordinates dual to (r, s) in T*X* and A(t, r, s, p, cr) is the pull- 
back of A(t, x, ^) under the change of coordinates (f, r, s, p, a) (f, x, ^) 



{t,r,s,dr(fi,ds<p), 
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We readily find that the equation D\fj{t, r) = implies 

dtifo{t, r) + A{t, r, s, dr(p, 4<p)|fl = 0, (4.4a) 

where dr(p\ii — dr'fio and dsip\R = J^Vid'Si. Analogously, from the real and imagi- 
nary parts of ds-D\ii{t,r) = we get 

^tv^ + a^, A + ^ dp^ Kdr^ Vt+Y^ da, Aip^j =0, ^ Axtj = 0, (4.4b) 

k j j 

and, finally, from the real and imaginary parts of dsiSjD\ji{t, r) we have 

dtipij + Ay -I- ^ ipiiBij + ^ *'Bu'Pi] + ^ VikCki^pij - ^ XikC'kixij = 0, 

/ / kA k,l 

9tXij + X! + X! + X! "f^kCkiXij + X! X'lkCkifij = 0, 

or, equivalently, 

9t<?;.y + Ay + ^ (l)aBij + 'Ba<l>ij + ^ (l)ikCki(jyij = 0, (4.5) 
/ ; kd 

where = Bji is the transpose matrix and the coefficients are 



A„(t,r) 



a^A ^ a^A a<y9j ^ d^K ^ a(/?j 



dsidsj ^ dsipk dvk ^ dr^ dpkdsj j-^ dph.dpk drh dvk ' 
„ , , d^A ^^A dip, ^ , ^ d^A 



9criasj ^ dcTidpk drk da^da^ 

In all the foregoing expressions, the arguments the the derivatives of A are evaluated 
at {t,r,0,dr(pa,J2Vidsi). 



Remark 4.3. Equation (^.5) constitutes a family of matrix Riccati equations and 
its unique solution exists in C°°([0,Ti]), < Ti < T. In addition, we require that 
(Xjj(0,r)) is symmetric and positive definite and it is not difficult to sfiow ttiat 
iXij) keeps symmetric and p ositi ve definite for t e [0,Ti]. Ttierefore, at least after 
redefining T, the solution of (4.5) exists in [0,T] with (xij) symmetric and positive 
definite. 

Remark 4.4. For the degenerate case R° — {x°} the coefficients take the form 
. . ^ 92A ^ , , d^A ^ , ^ d^A 



dsidsj ' ' dcTidsj ' ' daidaj ' 

which is the form known as beam tracing equations in the applied literature f4-V- 

In view of assumption^, the isotropic manifold A over the relatively open set O C 
Q, can be parametrized by coordinates (t, r), that is, points on An{(i, x, r, ^); (t, x) G 
O} are of the form {x{t,r),^{t,r)) with (t,r) e [0, T] x 7^°. 

Proposition 4.5. Let assumption ^ be verified. If we set (pQ{t,r) = i/j{x°{r)) 
and 'i(r,s)<P|i?X^, f) = C(^i^)i equations (4-.^) are identically satisfied. In addition, 
if (f>ij(t,r) solves the Riccati equation ( j^.jj j for t G [0,r] with initial condition 
'9siSj'0|if°(''); where ip is written in terms of the coordinates {r, s) near R° , then 
the complex phase ((.S) is a solution modulo 0{\s\^) of the Cauchy problem for 
the complex eikonal equation, that is, a) dt4> -\- \{t,x,dx4>) = 0(|sp), in O, and b) 
(f){0, r, s) - ipir, s) = 0{\s\^), in O". 
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Proof. The pull back of the vector field dt + d^X{t,x,dx(p) to the local coordi- 
nates (i,r, s) reads dt + d(^p ,j^A{t,r,s,drip,dsip), and the orbits of the latter field 
are just the orbits of the former field pulled back in the local coordinates. We 
know that the orbits on R oi dt + d^\{t,x{t,r),^{t,r)) in local coordinates are 
{(t, r, 0);t e [0,T],r = const.}, then dt + d(^p^^-)A{t,r,s,dr(p,dsf )\fi = dt and the 
partial derivatives 9pA, da A vanishes when evaluated for s = 0, (p,cr) = ^{t,r). 

On recalling that A and A are homogeneous functions of degree one in the dual 
variables ^ and (p, a), respectively, we can apply the Euler formula for which 

dtipo + Ht,r, s,drip,dsip)\s=o = dt^PQ + -^^^ + ^r~^3 = ^t'l^o = 0, 

^ dpk drk ^ da J 

which is solved by ipo{t,r) = claimed {ip\iio is real). Analogously, equa- 



tion ( 4.41: ) reduces to 

5tv?, + a,.A = 0, 

which is just the s-component of the Hamilton's equation d^/dt + d^X = expressed 
in the local coordinates, hence, it is identically satisfied. The only remaining equa- 



tion is the Riccati equation (4.5) which, indeed, is a nonlinear ordinary differential 
equation. If <pij(t,r) is a solution we have D{t,x,d(j)) — 0{\s\'^) and a) is proved. 
Analogously, b) follows from the Taylor expansion of the initial datum ip. □ 

The vice versa is also valid, that is, any approximate solution of the complex 
eikonal equation with the same initial datum should have the same reference man- 
ifold, determined by the characteristics, and the same Taylor polynomial around it 
up to the prescribed order; the proof requires the following lemma. 

Lemma 4.6. Let 5,tQ G 11+, / G (7(11+), and let us suppose that a non-negative 
function x G C"'^([0: ^o]; ) such that x(0) = and x{t) ^ ^ implies x'(^) ^ 
f{x{t))x{t)- Then, x{t) = for all t G [0,to]. 

Proof. Let = sup{ti G [0,to]; x{t) = when t G [0,<i]}; clearly t^ > 0. Since 
X G C^, if t* < to, there exists t'^ > t^ such that x(t) < 6 when t G [0,<'^]. We 
define k{t) = x{t) exp ( - J* f{x{t'))dt') and find k'{t) < in [0, i'J while fc(0) = 0. 
This implies x(i) = for t G [0,t'^] against the definition of i*, hence, the only 
possibility is = to- D 

Proposition 4.7. Let assumption Q he satisfied and (p ^ + ix & C'^(ri;C+) 
be a complex phase and let R = {{t,x) G il;x(/;,a;) = 0} be a submanifold, with 
(t, r, s) generic coordinate having the submanifold property for R. If (j) solves, mod- 
ulo 0(\s\^), the Cauchy problem for the complex eikonal equation with initial value 
ip satisfying assumption then, R — F{[0, T] x R°) and (f) has the Taylor expansion 
( l^.^ j with coefficients given in proposition |^.4 

Proof. Since Rr\X° — R°, we have dimi? — dimi?° -I- 1. According to the hypothe- 
ses, there are constants fci, fc2 such that |spfci < x ^ fc2|s|^ at least for \s\ small 
enough. Therefore, from the imaginary part of the complex eikonal equation near 
R, one finds 

dtx+ {dcX{t,x,dxip),dxX) <Cx^/^. 
We pick an integral fine t i-^ {t,x{t)) of the vector field dt -\- d^X{t,x,dxip) issuing 



from a point in R° and apply lemma 4.6 to x(i, x(t)) with the result that x(i, x{t)) 
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implying that the integral lines of V issuing from points in R° stay inside R. By 
applying the operator dx^ to the real part of the complex eikonal equation we have 

dtdx,(p + dx,X{t,x,dxf) + ^ Q^. dx- '^d^^^')- 

At the beginning of section ^, we have already argued that dxX\R = 0, hence, the 
latter equation implies 

Cvdxip\R + dx\{t, X, dxip)\R = 0, 

with Cv the Lie derivative along V . Therefore, along the integral lines of V we 
have the Hamilton's equations 

dx/dt — d^X(t, X, ^), d^/dt = —dxX(t, x, ^), 

where f = dx^p and d^/dt = ^^{Cvdxi<f)dxi. In view of assumption |4[ the solution 
x{t,x°), with initial position x° £ R°, exists for t G [0, T] and spans a manifold 
of dimension 1 + dim R° = dim R that should lie within R; we conclude that the 
whole reference manifold R is spanned by solutions of the Hamilton's equations 
and this justifies our construction. Finally, a straightforward application of the 
Taylor formula applied to D{t,x^d(f)) shows that the exact solution agrees with 
the asymptotic one ( |4.3| ) modulo 0{\sf') in the neighbourhood O <zVL where local 
coordinates [t, r, s) are defined. □ 



Proposition 4.5 allows us to achieve an asymptotic solution to the complex 
eikonal equation in a neighbourhood O D TZ where local coordinates {t, r, s) are 
defined. The relatively open set 7^ C i? is found in the form 7^ = F([0, T] x 7^°) 
with 7?.° C R" being the domain of definition for the local coordinates r. Therefore, 
one should address how to glue two solution based on coordinate defined in two 
different sets TZ", 7^2 1 with TZi fl 7^.2 7^ 0. This is straightforward since the normal 
vectors and, thus, the coordinates s, do not depend on the choice of coordinates 
r on R°. However, they behave non-trivially under a change of normal vectors e°. 

We conclude this section with an important lemma which means that a single 
complex phase (p can solve the complex eikonal equation corresponding only to a 
specific eigenvalue A in an open subset O C il. 

Lemma 4.8. // the complex phase (j> e C^{0) is a solution of dt4> + A(t, x, dx4>) = 
0(|sp) with X(t,x,^) being an eigenvalue of A(t,x,^), for any other eigenvalue 
Xi ^ X there are constants ci > and si > such that 

\dt(j) + Xi{t,x,dx(t>)\ > ci, for \s\ < si. 

Proof. Let us write 

\dt(t> + Xi{t,x,dx<p)\ = \{dt(t> + X{t,s,dx(l))) + {Xi{t,x,dx<i)) - A(t, s, d^^)) | 
> \dt(t) + X{t,s,dx(i))\ - \Xi{t,x,dx(t)) - X{t,s,dx4> 
We have \dt(f> + X{t,x,dx(l))\ < Ci\s\^, and, with Xi — {dxx)ij 

\Xi{t,x,dx(l>) ~ X{t,s,dx(t>)\ > \Xi{t,x,dx(l)) - X{t,s,dx(l))\ 

d^X 1^ d'^Xi 
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where \Xi{t,x,dx<j>) — X{t, s,dx(p)\ > C3 > in view of assumption ^ and 



dm 



'XiXj 



IE 



Pick < (Cs/Ca)!/^^ then, for \s\ < s,, 



dm 



'XiXj 



<C2\S\ 



\Xi{t,x,dx 
hence, 



X{t,s,dx 



lis 



d^x 



2 ^ 



■XiXj 



IE 



d^X, 



2 dm: 



'XiXj I 



|A/(t,x, (ia;^) - X{t,s,dx 



>C, = C^- C2SI > 0, 

1/31 



Analogously, if we pick s; < min{s,, (C,/Ci)^''^}, we have Cisf < C, and, thus, 



5t0 + A(t, s,dx 



< C*, and for \s\ < s;, 
|5t</) + Ai(i,a;,d^(^)| > - Cisf = q > 0. 



□ 



Remark 4.9. A shorter proof of lemma ^.S is obtained on noting that, if the claim 
is not true, for every C > there is A; such that |A — Ajj < C at least in a point 
on A against assumption The extended proof given above is somewhat more 
informative. 

5. The Amplitudes 

Let (f> be an approximate solution of the complex eikonal equation addressed in 
section ^ with A a generic eigenvalue of A corresponding to the projector tt; in this 
section we shall address equations (2.7) and (2.8) with 0^ = and q — 2; upon 
writing a|[''' — oq and a'j}"' = ai for simplicity, we shall study the equations 



aLo{t,x,d(t))ao{t,x) = 0{\s\^), 
aLo{t,x,d4>)ai{t,x) + L{t, x, d)ao{t, x) = 0{\s\), 



(5.1) 
(5.2) 



in a neighbourhood of i? = {{t, x) e fl; Im(/)(i, x) = x{t, x) = 0} where coordinates 
{t, r, s) are defined; we recall that, according to proposition (4.7), R is the flow out 
oi R" — R n {t = 0} along the characteristics of Lq. 

We have ag = (/ — vr)ao + ttoo, with tt the extended projector (with n — 2), and 

aL,X^,x,d(i>)TTaQ = i{dt(l) + X{t,x,dx(i>))T^aQ + 0(\sf) = 0{\s\^), 



in view of corollary 3.4. Therefore, 

aLo{t,x,d(j>)ao = crL„{t,x,d<j)){I - 7r)ao + 0{\s\^) 

= i{dtcl) + Xi{t,x,dx(l)))niao + 0{\s\^), 



and, on applying lemma 4^, we have that equation (5T) is equivalent to 

(/-i^)ao = 0(|s|3). (5.3) 



Analogously, for the inhomogeneous equation ( |5.2| ), we find that the component 
TTOi is arbitrary since crL(,7rai = 0(|sp), so we can set nai — and ai — {I — 7r)ai. 



16 



OMAR MAJ 



Then, by applying the extended projectors tt and / — tt to the inhomogeneous 
equation we get the necessary condition 

n{Lo{t, X, d)ao + B{t, x)ao) = 0{\s\), (5.4) 

since iruLoai — 0{\s\^), and the algebraic equation 

(/ - ^)ai„ (t, X, d<j)) (/ - i)ai + (/ - {Lo{t, x, d)ao + B{t, x)ao) =0{\s\). (5.5) 

Equations ( ^.3|) -( ^r5| ) should be satisfied in the neighbourhood O where coordinates 

(t, r, s) are defined. 

On considering first equation ( |5.3| ), we write 

Tr{t,x,dx(f)) = 7r(i,r) + ^7rj(t, r)si + ^^7rij{t,r)siSj + 0(|sp), 

where 7r{t,r) = 7:{t, x, dx(p{t, x))\r, whereas TTi{t,r) = 9s,7r|fl and nij = Os^s^t^Ir- 
The Taylor expansion for the amplitude show that equation (5.3) determines the 
first three coefficients of the expansion only. Clearly, the space of second-degree 
polynomials in s with coefficient s in C°°{R) constitutes a C°°(i?)-module. We 
now show that the solutions of (5_^) are one-to-one to functions in C°°(i?;C^) 
satisfying the appropriate polarization condition. It is worth noting that 7r(t, r) is 
the projector TT{t,x,^) restricted to the isotropic manifold A, locally parametrized 
by it,r). 

Proposition 5.1. Afunctional £ C°°(ri;C^) solves (5.S) if and only if its restric- 
tion to a neighbourhood O of R amounts to M{t,r,s)a{t,r) where a £ C°°(i?;C^) 
satisfies the polarization condition {I — 7r(<, r))a{t, r) ~ and 

Af = / + V Ah{t, r)s, + ^ V My(t, r)s^Sj, 



Mi,M,j e C°°(i?; End(C^)), TT{t,r)M,{t,r)n{t,r) = Tr{t,r)lVUj{t,r)Tr{t,r) = 0. 



First, it is useful to prove the following lemma. 

Lemma 5.2. With the notations given above we have Tr{t,r)Tri{t,r)Tr{t,r) — and 

TT{t,r){TTi{t,r)Tij[t,r) + ■nj{t,r)Tii(t,r) + ■Kij{t,r))TT{t,r) = 0. 

Proof Since _ ^ ^ 0(|sp), i.e., ds,i^^ - n)\R = and dsAA^^ " ^)\r = 0, 
the claim follows on performing explicitly the derivatives and multiplying by 7r(t, r) 
both on the right and on the left. □ 



Proof of proposition 5. 1 . Let a, , be the coefficients of the expansion of oq near 
R modulo 0(|sp). Then, we have 

(/ - 7r)ao = (/ - 7r(t, r))a + ^ [(-^ - 7r(i, r))ai - ■Ki{t,r)a]si 

+ [(^ - 7r(i, r))ay - TT^{t,r)aj - 'Kj{t,r)a^ - TTij{t,r)a]siSj -f 0(|sp). 

We see that ag solves equation ( |5.3| ) if and only if 

(/ - TT{t, r))a = 0, (/ - 7r(t, r))ai = 7ri(i, r)a, 

(I — 7r(t, r))aij = 'Ki{t, r)aj -f 7rj(t, r)ai + T^ij{t, r)a. 

The first equation is identically satisfied in virtue of the hypotheses. Lemma 
5.2 provides the solvability conditions for the remaining equations the solution of 
which is thus readily found with the result that Mi{t,r) = Tri{t,r) and Mij{t,r) — 
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TTi{t,r)Trj{t,r) + ■nj{t^r)'Ki[t^r) + Trij{t,r). Finally, again by means of lemma 5.2 



one can see that 7r(i, r)Mi {t, r)Ti{t, r) — 7r{t, r)Mij (t, r)Ti{t, r) —0. □ 
There is a natural smooth extension of Ala to the whole domain f2 given by 

where a(t, x) is any smooth extension of a{t, r) to a compact neighbourhood of R, 
Xi{t,x) = dx{t,x)/dxi and the derivatives of tt are to be evaluated at {t,x,dx(p). 
This follows by direct substit utio n into (5_^) and by using the following identities 



that can be proved as lemma 5.2 



Lemma 5.3. If tt{x) is a projector-valued -function of x G R'' , then TrSiTTTr = 0, 
and TT^diirdjir + djirdiir + dijTr)Tr — 0, where di — d/dxi and dij = /dxidxj . 



Let us now consider ( p.4[ ) with oq given above. First we note that, 
9(00 = dtina) + 0{\s\), d^^aQ = 3^^(1:0) + i'^d^.Trira + 0(|s|), 

hence, 

Lo{t,x, d)ao = Lo{t,x, 9)™ + dx dx ^'%^^)-+ '^d'^D' 

and equation reads 

'KLo{t,x, 8)110, + 2 ( Z_/ Qrj. Qrj. ^(^»%^ + ^j%7r)7rja + TrSTra = 0(|s|). 

The differential term in this equation can be further simplified on writing 

nLona = TiCva + 7r(Lo7r)a, 

where where Lv is the Lie derivative along the vector field V = dt-\- d^X{t, a;, dx^p) 
and we have used the identity 

dA fdiAir) dn \ / dX \ ^X 



TT^feTT = TT „ ^ TT = TT 



.OTTX (OA ^ , .^OTT \ OA 



96 V 96 

As for the second term, the symmetric part of Aid^^i: can be computed as follow. 
From one hand, we have 

9^(A7r) (dAdiT a4 97r\ d^i: 

where A = ^ A^^^, thus, Ai — dA/d^i and A/ dS^id^j = 0. On the other hand, 
d^jAn) _ d^jXir) _ d^X /^AStt SASttx d^n 

so that 



Finally, in virtue of lemma pTsl wc find 



d'^X 
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Since V is tangent to R we find that equation (5^) is satisfied if and only if a{t, r) G 
C°°(i?;C^) solves the transport equation 

J (/-7r)a = 0, 

I -kCvcl + i:{Lqi: + B)a + ig^a = 0, ^ 
where all the coefficient should be evaluated on R and 



2 dxidxj d^id^j 

^^3 



introduces a phase-shift in the amplitude due to the wave field localization. Apart 
from such a phase-shift effect this is exactly the geometric optics transport equation 
evaluated on R. We also note that, since 9ao|fl is independent on the choice of the 
smooth extension oq, the transport equation is also independent on such choice. 



The transport equation (5_^) allows us to obtain a solution for a by integrating 



along the integral lines of the vector field V , that are parametrized by R° = Rn{t ~ 



0}, then, we get a solution f or a p given in terms of a by proposition 5.1 . 

Finally, we have to solve ( |5.5| ) which is purely algebraic and one readily see that, 
with ai = (/ — 7r)ai, it is equivalent to 

[(/ - 7f)crLoai -I- (/ - 7r)(Loflo + Baa)] = 0. 

Since 

{I-n)aLo^ ^ i{dt(f) + \i{t,x,dx(l)))T:i{t,x,d^(j)), 
i;A,#A 

lemma ensure that such an operator is invertible on a neighbourhood of R and 
denoting by Q{t,r) its inverse evaluated on R we can set ai(t,x) equal to any 
smooth extension of 

- Qit, r) [{I - n)iLoao + Bao)] (5.7) 
to a compact neighbourhood of R. 

6. Construction of the Complex Geometric Optics Solution 

Now we can construct the complex geometric optics solution of the Cauchy 
problem for (2.1) provided that assumptions hold true together with condition 



First, we have to solve the Hamilton's equations and obtain the isotropic mani- 
folds Ai, . . . A,„ along with the corresponding projections . . . , Rm that give the 
reference manifold R — [J^Rfi] the existence of such geometric objects is ensured 
by assumption ^ in section |^. 

Then, we make use of the construction of section |4| in order to obtain the coordi- 
nates patches Of^j on the basis of a finite covering {T?.^,^ }£ of i?° ; that exists since, 
in view of assumption ^, i?° is a closed subset of the compact set X°, hence, it is 
compact. In each neighbourhood 0^_f we have coordinates {t, r, s) G [0, T]xOr'xOs- 

Next, we obtain the approximate solutions 0^ to the Cauchy problems, 

dt<P^, + Xi{^){t,x,d^(t)f,) =0, 0p|t=o(a;) ^ipf,{x), 
with Z(/i) being given in condition |l|. In the statement of proposition |2.5| we have 



assumed djxfj.\R°{x) > and, according to remark 4.3, it is X/j(^j?':s) > 
for some constants > 0: this allows us to apply lemma with q = 2. The 
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phases (f>^ are defined in = IJ^ O^^^ up to a remainder of 0{\s\^), therefore, they 
are better understood as a representative of an equivalence class for the following 
equivalence relation in C°°(il; C™), 

every two <j), (f) E C°^'{il; C™) are equivalent if and only if each component 
(pfi and 4>'^ have the same second-degree Taylor polynomial in the variable s 
near the submanifold i?^. 

We pick a representative 4>{t,x) such that Im^^ > outside IJ^ O^^f. Given 
such a representative, we can apply the results of section ^ in order to obtain a 



solution al,"'' and aj^^-* of ( ^.7| ) and (2^) with q — 2. Such solutions are defined 



in the neighbourhood modulo a remainders 0(|s| ) and 0(|s|), respectively; in 
analogy to the complex phases, this defines an equivalence relation in C°°{^1; C^) 
of functions with the same Taylor polynomial near i?^. 

Finally, we can construct the (lowest order) complex geometric optics solution 



(2.5) which is therefore defined modulo the foregoing equivalence relations. 



Proof of proposition 2.5. The final part of the argument relies on the Maslov's esti- 



mates proved in lemmas 2.2 and 2.4 as well as on a partition of unity subordinated 



to the finite covering of R° with open sets 7?.° ^ C R° . 

For each open set 7^° ^ let us consider the neighbourhood O^^ C fl where coor- 
dinates (i, r, s) are defined. Assertion (iv) of assumption ^ allows us take O^^^ so 
small that Of^j n O^^f- = if ^ ^ The open sets 0° ^ = Of^j Ci X° satisfies the 
hypotheses of the partition of unity at R° Theorem 1.4.5]. Therefore, one can 
find functions ^ G C(5"(0° ^) such that a;° ^ > and £ w° ^ < 1 with equality 
in a neighbourhood of R° . 

By using coordinates (i,r, s), we define the functions ujf^{t,r,.s) — Lu°{r,s) and 
we see that ujf^j € C°°{0^^i) and, for every t € [0, T], 

LUf^Atr) e C^iOf^jnX') and 
ujf_i^i{t, •) = 1 in a neighbourhood of i?^ = R^ n X*. 



The sum is over i only as, in a neighbourhood of i?^, we have ujiyj = Q li v ^ 

Proof of a). Let us write 1 = (l — ^^.f) + Y^i i so that (1 — Yli i) is 
supported away from i?° whereas a;° ^ is supported in 0°^. More specifically we 
have. 



and, in virtue of lemma the last term is < C^^^s^ for every A: S IN, with the 
sup being computed over S = supp(l — '^gOj"^ g) C X° . As for the first term we 
exploit the Taylor's formula for il)^ and, on recalling that ^^it^o amounts just to 
the second degree Taylor polynomial of i/;^, we get 
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where the remainder p^(r, s) has been obtained from the fundamental theorem of 
calculus, 

|q|— 3 

and the integral is uniformly bounded by 1. Therefore, we have 

1/,)° (h p''V'n/e _ „e \\ < \(,,° (h — n^^^ \ pi4> ^\t=o / \ 



The first term is of ©(e^) in virtue of lemma 2.4 applied at t = and with the 
compact set K = supp(ijj° ^), as and a^'JI^g are the same when evaluated on 
The second term is 0(e) uniformly on whereas the third term can be estimated 

by 

|w°,££'m^"^"'*^°^1 ^ - I] |s"e*'^"'*="^1sup|w°^<,V;.,a| <e^ const., 

^ |a|=3 

again, uniformly in X° , where we have used the fact that 

|sae#,,|t=o/e| < ^||^ag-™^|^ i; = s/£5, |a|=3, 

and this proves assertion a). 

Proof of b). We write Lf^ = (1 — a;^/)if;^ + w^^^Li;^, where L = Lq + 
and we have, cf. equation 

with 

i 

uniformly in Q and for every /c G IN in view of lemma ^.2| . Moreover, 



We know that equations (2.7) and ( |2.8| ) hold so that the second degree Tay- 
lor polynomial in s of lo ^^a Lgolf^^ and the zero-degree Taylor polynomial in s of 
io ^(a L^a'j}'' + L^af^^ + SaJ^') vanishes. Hence, we can apply lemma 2.4 on the 
compact set supp(w^,£) C 51 so that the first two terms are 0(e2 ) uniformly on the 
whole domain il. Finally, the last term is 0(e). □ 



Proof of corollary \2.(\ . For every fixed e £ R+, we know that there exists a unique 
solution G C°°(ri), which satisfies the energy inequality, 

\\u'mLHxn<e'''\\u'mLHx^), 

where the constant K depends only on L°°-norms of J2 dxi^i{t, x) and Bit, x). In 
particular, K does not depends on e and {u^(0);e G (0,eo]} is a bounded family 
in L^(X°), thus, G (0,£o]}, with u'^{t,x) extended by zero outside Q, is a 

bounded family in C{[0,T], L^(R'')). 
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Let v"^ be any representative of the equivalence class in proposition 2.5 



m a 

common domain of determinacy il. Then, the difference = u'^ — satisfies the 
problem 

r L{t, X, d)w^t, x) = r{t, x) e C°-{Tl- C^), 
\w\t=o{x) ^ g^{x), 
where, in view of proposition ^.5|, 



\\9 



=sup|/(x)| <Ce^ liriUoo(o) =sup|r(t,a;)| <C'£i 



In addition, for < t < T, we have the inhomogeneous version of the energy 
inequality 

< e*^ll/llL^(xo) + re^*-*')^Iir(OIlL.(x*')^^'' (6-2) 
Jo 

with the same constant K as before. Then, we find 

\\w%t)\\L2(xr) < a(t)|l5llL-(f2) +fo(i)||/1lL-(o) < cit)eK 
and c G C([0,r];]R+). □ 
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